A quantum decay model with exact explicit analytical solution 
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A simple decay model is introduced. The model comprises of a point potential well, which 
experiences an abrupt change. Due to the temporal variation the initial quantum state can either 
escape from the well or stay localized as a new bound state. The model allows for an exact analytical 
solution while having the necessary features of a decay process. The results show that the decay is 
never exponential, as classical dynamics predicts. Moreover, at short times the decay has a fractional 
power law, which differs from perturbation quantum methods predictions. 
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Despite the fact that almost all the decay processes in 
nature observed of having an exponential decay law, it is 
well known that quantum mechanics predicts deviation 
from this law [H, y, Q • 

In fact, it has been proven by Khalfin 4] that the long 
time behavior of the non-decay probability cannot be ex- 
ponential, and in practical cases it has a negative power 
decay law. Moreover, at short times most quantum sys- 
tems obey t'^ dependence. This kind of dependence is 
usually attributed to a reversible process, which contra- 
dicts an irreversible decay to the continuum. 

Recently, it became technologically feasible to measure 
this deviation from exponential law, and indeed it was 
demonstrated experimentally [g, S, 0] ■ 

The scenario of an irreversible decay of a confined 
bound state to the continuum is ubiquitous in the phys- 
ical world. Beta decay is such a realization. Hence, re- 
solving this controversy is required to the understand- 
ing of these basic processes. Theoretically, this problem 
was confronted by applying an abrupt perturbation on 
the initial confined state. However, this perturbation ap- 
proach merely emphasizes the controversy except for the 
intermediate times where an approximately exponential 
law seems to appear [8|. At short times the reminiscent 
of a reversible t^ law is still dominant. 

The recent technological developments, which allow 
trapping cold atoms in very small traps [3, [l^] , also al- 
low to release them almost instantaneously (since the 
trapping is done by laser beams). As a consequence, 
the temporal dynamics and decay of quantum particles 
at the presence of an abruptly changing potential can 
also be investigated in the laboratory. In this paper we 
investigate a simple quantum mechanical model, which 
can emulate realistic decay scenarios. The initial state 
is a bound eigenstate of a localized potential well, and 
the model allows investigating the state dynamics due 
to abrupt change in this potential well. To simplify 
the model we use a delta-function potential well. It is 
well known that this kind of a potential can emulate any 
barrier/well whose de-Broglie wavelength is considerably 



longer than the potential physical dimensions |ll| . In 
other words, for most practical purposes it can replace 
any point potential. On the other hand, a point poten- 
tial brings in singularity to the system and as a result 
a fractional power law emerges. However, as was men- 
tioned elsewhere [l2| some reminders to this behavior can 
be traced even in the analytic potential case. 

The main strength of this model is that it has an ex- 
act analytical solution and no approximations are taken. 
We show that the dynamics of this model does not look 
exponential at any time. Moreover, at short times (as 
well as at very long times) the dynamics is governed by 
a fractional power law instead of an integral power law. 
To the best of our knowledge, this is the only quantum 
decay model, where the final state can be either localized 
or extended, and has an exact analytical solution. 

We use a delta function to simulate the attractive po- 
tential well 



V{x,t) = 



-2ad{x) 
-2\S (x) 
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That is, the potential can be written as an initial sta- 
tionary potential Vq (x) = —2aS (x) and a perturbation 
part AV {x, t) = 2 (a — A) u (i) S {x) [u (x) is the heav- 
iside function u (x) = {0 for a; < 0; 1 for x > 0}]. The 
Schrodinger equation reads 



i^ = -^ + Vo{x)i; + AV{x,t)^ (2) 

hereinafter we adopt the units h ^ 1 (Planck constant) 
and 2to = 1 (particle mass). It should be stressed that 
any shallow potential well {Ax >> (C/q) with width 
Ax and depth — C/q) can be replaced, for most practi- 
cal purposes, by a delta function potential Vmiuai (x) = 
—2a6{x) with the prefactor 2a = UqAx, since the two 
scatterers have a very similar scattering and a single 
bound state ( for a positive potential see, for example 

iUi])- „ 

By applying the same analysis and logic of [13| it can 
be shown that when the well has finite width (instead of 
a delta function), say Ax, then all the derivations that 
follows are valid provided t >> Ax^. Therefore, this 
behavior can be traced for every shallow potential well. 



i.e., where its eigen-boundstate energy Eq obeys \Eo\ « 

Next, for simplicity, we renormalize space and time to 
the dimensionless variables Xnew = ctx, tnew — to^and 
choose the normalized parameter /i = A/a. With these 
dimensionless parameters, the potential can be written 



T/(x,t) 



-26 (x) for t < 
-2^5 (x) for t > 



Therefore, if the initial state ilji{x,t) is the bound 
eigenstate ipBi {x, t) of the unperturbed well, then 
^pi (x, t) = ipBi {x, t) = exp (— \x\ + it). After the abrupt 
potential change, the only localized state is, of course 

^Bf {x, t) — fi^^^^ exp [~fi \x\ + ifi^t) . (3) 

This model has an exact solution without the need 



for any simplifying approximations. The wavefunction 
at any instant can be calculated by the integral expres- 

oo 

sionip{x,t) — J K {x,x';t)ipi{x' ,t — 0)dx' where the 

— oo 

Kernel of the integral (the Green function) is [13] 

K (x, x'; t) = Kfree {x, x' ] t) + 



— exp [— ^ (|a;| + \x'\ — i^t)] erfc 



\x'\ -z2/it 



2\/it 



(where erfc is the complementary error function [141] ) 
and the free space Kernel is 

. ix — x') 



^ free [Xj X , t) 



■- exp 



After some tedious, albeit straightforward calculations, 
the solution for the initial wave function (ipBi {x, t)) can 
be written 



i^{x,t) 
1 



2'' 



■1=^1 erfc 



2Vii_ 



jl^'erfc 



2V^ 



r 



(4) 
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It should be stressed that Eq|4]is the exact solution with- 
out any approximations. 

Note that for /x = 1 there is no change in the poten- 
tial, and therefore Eq. |4] degenerates to the eigenfunction 
ip {x,t) = ipBi {x,t) = exp (— |x| -I- it), which means that 
the wavefunction remains in its initial state. 



On the other hand, when /i = 0, after the transition 
(the abrupt change) the potential vanishes and the wave- 
function propagates freely in space 



^p{x,t) 
1 , 



(5) 
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FIG. 1: The distribution of the probability density in space 
for four different times: t = 0.07,0.2,0.7 and t = 100. The 
green dashed line represents the initial states \ip{x,0)\ , the 
red dash-dotted line is the final bound state \^{x,oo)\ , and 
the solid blue line is \^p {x,t)\ . 



This is a relatively simple but important case (due to 
its generic nature), so we would like to elaborate on its 
dynamics. 

In this case If/) (x, 0)| and {ip {x, t)\ are similar (except 
for the oscillations) only till a certain x, beyond which 
the wavefunction decays like |'0(x,t)] ~ x~^ since for 
x^/t » 1 

; / ,\ r-^ 1 exp(ix^/4t) 

This result is consistent with the prediction of ref . [I4] 
that the wavefunction at very short times and long dis- 
tances, i.e., x^/t >> 1, is 



V'(x,t)~[^'(o+,o)-v/(o-,o)] 



■ exp 



m 



(the initial exponentially small value of the wavefunc- 
tion at 2: ^- 00 was ignored). When < /x < 1 the 
dynamics is more intricate since the final states can be 
either extended (as in the /i = case) or localized (Eq.2). 
The plot of the probability density \tp {x,t)\ as a func- 
tion of X, as depicted in Fig.l illustrates this point. 

When the perturbation is turned on the initially lo- 
calized particle's energy is modified and the particle can 



remain localized at a different energy, i.e., Eq = —fj? (in- 
stead of the initial one E^ — —1) but it can also escape 
to the continuum. 

At short times, t << x^ the wavefunction can be ap- 
proximated by 



4> {x, t) = e 



^ „it-\x\ 



(itf^ exp iiu^/U) 
4(1-m)^^ ' ^ , ^ (6) 



For X — > oo the wavefunction's pertrubative term decays 
like x^"^ . In EqlHlwe see that the short time behavior have 
a fractional power law and deviates from the reversible 
t^ dependence Q. 



In the long time regime, i.e., t 
can be approximated 



oo, the wavefunction 



V' (a;, t) ^ {^1-^ - 1) W - |a;| r^l^ exp {ix^/U) + 



4/x 

1 + /X 



exp (z/i"i — /i \x\) . 




FIG. 2: r/ie temporal revolution of the probability density for 
fj, = 3 (solid line). The dashed line stands for its final value 
P{t = oo). 



We observe two different dynamics regimes. At very 
large distances from the origin (but still i >> a;) the first 
term, which is related to the propagating-waves rules, 
while at short distances the first term is merely an os- 
cillating correction to the second one, which is related 
to the final localized state. As i ^ c» the wavefunction 
converges to the final bound state ^ with extra factor 

Usually, the measured quantity, which quantify 
the decay rate, is the survival amplitude A (t) — 

oo 

J ip* {x,0) ip {x,t) dx] and the survival probability 

— C30 

P (t) = \A {t)\ is the probability to remain in the initial 
state (i.e., the non-decay probability); similarly, 1 — P (t) 
is the probability to escape to infinity, i.e., to decay. 

A (t) can be calculated exactly and straightforwardly 
(albeit with tedious calculations). For the initial state 
{■>pBi {x,t)) we find: 



J 



A{t) 



{l + f^Y 



I erfc (Vit) [l + /? - 2it (l - ji"^)] + 2 (l - n"^) J-e 



"**-H2^e**(^'-i) 



2 — crfc 



I yitiJ. 



(7) 



For completeness we add the special ^ — case: 



A (t) = erfc ( Vit ) [1 - 2it] + 2y -e 



(8) 



In Fig. 2 the dynamics of the survival probability is 
presented for fi = 3. It is clear from the figure that the 
probability decays eventually irreversibly to a constant 
value. However, it never decays exponentially. 

At long time the survival amplitude goes like 



A{t^oo) 



4/ie 



it(/j^-l) 



(i+Mr 



4^^1/2 (i)3/2 



^ — itfi — «3'7r/4 



and the non-decay probability can be approximated 



P{t) = \A{t ^oo) 



(1 + m)' 



1 + ^^ ^ ;,, cos (t^l^ + 37r/4) 



It oscillates with angular frequency iJ? with varying 
amplitude that decays like t~^/^ and converges to the 

value PW-(^. 

It should be noted that this final probability is smaller 
than 1 for either /^ > 1 or ^ < 1. This result obvi- 
ously contradicts the classical intuition that the non- 
decay probability decreases only when the well is raised. 

Despite the irreversible nature of the process, it has no 
similarity to the well-known exponential decay. 

At short times i << 1 this expression can be expanded 
by fractional powers series to 



Therefore the non-decay probability goes like 



A(t « 1) - l + 2zt(^- 1) 



;(_l)3/4(^_l)2^3/2 



(m-1) Mt' 
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FIG. 3; The escape probability 1 — P (t) as a function of time 
(the solid black line) for /i = 3. The dashed red line stands 
for the short time approximation (Eq\§j). 



l-P{t«\) = l-\A{t« 1)1' 



(9) 

which means that the leading term in the escape prob- 
ability has a fractional power dependence on time (see 
Fig.3). 

This behavior resembles l.3| where the dynamics of an 
ad-hoc potential spectrum was investigated; however, one 
of the advantages of the model present here is its physical 
realization. To the best of our knowledge, this model is 
the only decay model, which allow for an exact explicit 
analytical solution. 

The /i ^ 1 regime calls for comparison with perturba- 
tion methods, which lead to the Fermi Golden Rule. In 
the latter case the short time regime goes like t^ instead 
of i3/2 of EqU 

As was said at the beginning of the paper, it should 
be stressed that even if the well had a finite width Aa; 
(instead of a delta function one) , then the fractional be- 
havior, which appears at Eqsl7][9] would still be traced for 
every shallow potential well (i.e., \Eq\ ~ {Uq/S.x)'^ / A « 
Ax~^), and it is not merely a mathematical anomaly. 

To summarize, the dynamics of a perturbed delta func- 
tion potential well was investigated. Although this sce- 
nario can model a realistic case (such as a particle decay 
from a point potential trap), it has an exact analytical so- 
lution. Not only does this model behave differently than 
the well-known exponential decay law as classical decay 
laws predict, but it does not even have an integral power 
law at short times as quantum processes suggest. In fact, 
the dynamics is more intricate and has a fractional power 
law at short times. 

We believe that the analyticity of the solution of this 
model along with its experimental feasibility can be used 
to shed light on the generic decay process from both prac- 
tical as well as theoretical perspectives. 



[1] L. Fonda, G.C. Ghirardi, and A. Rimini, Rep. Prog. Phys. 

41, 587 (1978) 
[2] P.T. Greenland, Nature 335, 298 (1988) [9] 

[3] C.B. Chiu, E.G. G. Sudarshan, and B. Misra, Phys. Rev. 

D 16, 520-529 (1977) [10] 

[4] L.A. Khalfin, Zh. Eksp. Tear. Fiz. 33 1371 [JETP 6, 

1053-1063 (1958)] 
[5] V.J. Menon and A.V. Lagu, Phys. Rev. Lett. 51, 1407 [11] 

(1983) [12] 

[6] S.R. Wilkinson, G.F. Bharucha, M.G. Fischer, K. W. 

Madison, P.R. Morrow, Q. Niu, B. Sundaram and M.G. [13] 

Raizen, Nature, 387, 575-577 (1997) [14] 

[7] P.T. Greenland, Nature 387, 548-549 (1997) 
[8] See, for example, C. Gohen-Tannoudji, B. Din and F. 



Laloe, Quantum Mechanics Vol. II (Hermann, Paris 
1977) 

P. Szriftgiser, D. Guery-Odelin, M. Arndt, and J. Dal- 
ibard, Phys. Rev. Lett. 77, 4 (1996) 

G. Fort, L. Fallani, V. Guarrera, J. E. Lye, M. Mod- 
ugno, D.S. Wiersma, and M. Inguscio, Phys. Rev. Lett. 
95 170410 (2005). 

E. Granot, Phys. Rev. B 71, 035407 (2005) 
E. Granot and A. Marchewka, Europhys. Lett. 72, 341 
(2005). 

M.Kleber Phy.Rep. 236,331 (1994) 

M. Abramowitz and A. Stegun, Handbook of Mathemat- 
ical Functions , Dover publication (New- York 1965) 



